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Abstract 

In a recent study [T, an extension of the M/M/l queueing process with the excluded- volume effect as 
in the totally asymmetric simple exclusion process (TASEP) was introduced. In this paper, we consider its 
discrete-time version. The update scheme we take is the parallel one. A stationary-state solution is obtained 
in a slightly arranged matrix product form of the discrete-time open TASEP with the parallel update. 
We find the phase diagram for the existence of the stationary state. The critical line which separates the 
parameter space into regions with and without the stationary state can be written in terms of the stationary 
current of the open TASEP. We calculate the average length of the system and the average number of 
particles. 

1 Introduction 

The queueing theory has been considerably studied since Erlang started designing telephone exchanging system 
in 1909 [S\. In [9], he developed the theory of the call-loss system, which significantly contributed to the progress 
of telephones and electric communication systems. After his study, Kendall presented his first paper [15]. Since 
then, the study on the queueing theory including Kendall's notation [TB], Burke's theorem [3], Jackson networks 
[13], and Little's theorem [TH] has been accelerated. Nowadays, it is applied to study on social systems such as 
the Internet [20] , resource management systems [5] , vehicular traffic systems [T2] and pedestrian traffic systems 

[TO 

Let us review here the discrete-time M/M/l queueing process, which will be extended in the next section. 
M/M/l means that the queue has a Markovian entry of customers and a Markovian exit at 1 server. Each 
particle (customer) enters the system with a probability a and receives service with a probability 0, see figure 
[T] The system is characterized only by the number N E Z>o of particles. The probability P{N; t) of finding N 
particles at time t is governed by the following master equation: 

P(0; t + 1) =(1 - a)P(0; t) + a0P(O; t) + (1 - a)0P(l; t), (1) 
P(N; t + 1) =a(l - 0)P(N - 1; t) + (1 - a)(l - 0)P(N; t) 

+ a0P(N;t) + (l-a)0P(N + l;t), ^ ' 

for N e N. When a < j3, a unique stationary-state solution to the master equation can be easily obtained as 
the following geometric distribution: 

1 > 0(1 -a) \(l-a)0j {6) 
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Figure 1: M/M/l queueing process. 
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The average number of particles can be easily calculated as 
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On the other hand, the asymmetric simple exclusion process (ASEP) in one dimension is one of typical 
interacting particle systems and admits exact analyses of non-equilibrium properties [181 124] . It has been 
vigorously studied in the recent two decades. The stationary state of the totally ASEP with open boundaries 
(open TASEP) was found in a matrix product form in [7 . Since then, matrix product stationary states of 
various generalizations of the ASEP containing discrete-time versions of the open TASEP have been found [3]. 
The TASEP is one of the basic models of the vehicular traffic [B] . 

Although both the M/M/l queueing process and the open TASEP are Markov processes with particle entry 
and exit, they are different in the following two points. First, the number of particles in the M/M/l queueing 
process does not have an upper limi10, whereas that in the open TASEP cannot be greater than the fixed 
number of sites. In other words, the state space of the M/M/l queueing process is infinite and that of the 
open TASEP is finite. It should be noticed that the open TASEP is a call-loss system, which has been seldom 
discussed. Second, the TASEP has a spatial structure, and particles are affected by their excluded volume 
(hard-core repulsion). By contrast, the M/M/l queueing process has no spatial structure, and only the number 
of particles characterizes the system. 

An extension of the M/M/l queueing process with the excluded-volume effect (in other words, the TASEP 
with a new boundary condition) was introduced in [T]. In this model, particles enter the system at the left site 
next to the leftmost occupied site and exit at the rightmost site of the chain, see Fig. [5] When we (pedestrians) 
make a queue, we usually proceed if there is a space in front of us (excluded-volume effect), which is a motivation 
to treat this model. Let us call the model exclusive queueing process. 

In pP , a stationary state of the exclusive queueing process was given by a slightly arranged matrix product 
form for the open TASEP, where the probability of finding each configuration is proportional to 



The row vector, the matrices and the column vector are independent of the entry rate a, whereas c is independent 
of the exit rate /?. The exponent L denotes the position of the leftmost particle. 

In this paper, we consider a discrete-time version of the exclusive queueing process with the entry, hopping 
and exit rates (a,p and /3, respectively) replaced by probabilities within one time step. Although there are some 
choices of an update scheme, we take the fully-parallel-update rule@ The special case where the bulk hopping 
is deterministic (i.e., the hopping probability is 1) has been analyzed in |25j . 

To find a stationary-state solution is one of basic problems. The idea to find it in this paper is reducing 
the balance equation for the exclusive queueing process to that for the open TASEP. Then, a stationary-state 
solution will be obtained in a similar form to that for the continuous-time model as ([5]). We calculate the critical 
line which separates the parameter space into the regions with and without the stationary state by evaluating 
the convergence of the normalization constant. The critical line will be written in terms of the stationary current 
of the open TASEP. We calculate the average length of the system (the position of the leftmost particle) and 
the average number of particles in the system. We also calculate the average number of particles in the queued 

This paper is organized as follows. In Sect. [21 we define the exclusive queueing process with discrete time. 
In Sect. [31 we obtain a stationary state of the model. In Sect. [H we identify the region where the stationary 
state exists in the parameter space, and calculate the average length of the system and the average number 
of particles. In Sect. we treat the model in the case where the bulk hopping probability is 1. Section [6] is 
devoted to the conclusion of this paper. Detailed calculations are done in Appendices. 



We introduce an extension of the discrete-time M/M/l queueing process with the excluded-volume effect on the 
semi- infinite chain (exclusive queueing process), see Fig. [21 where each site is labeled by a natural number from 

1 Queueing processes with an upper limit of the number of particles have been studied. 1 141 

2 The discrete-time open TASEPs with the random-, forward-sequential-, backward-sequential-, sub-lattice-parallel- and fully- 
parallel-update schemes have been studied, see [3] for review. It is reported that the movement of pedestrians in one dimension is 
well represented by the parallel- update dynamics 1221 . 

3 We define the number of particles in the queue by counting particles except for one at the rightmost site. 



c L x (row vectorj matrix product |column vector). 



(5) 



2 Model 
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Figure 2: Exclusive queueing process. 



right to left. Each particle enters the chain at the left site next to the leftmost occupied site with a probability 
a, hops to its right nearest neighbor site with a probability p, if it is empty, and exits at the right end of the 
chain with a probability f3. If there is no particle on the chain, a particle enters at site 1. These transitions 
occur simultaneously within one time step. In other words, we take the fully-parallel-update scheme. (We call 
it simply parallel update hereafter.) The model is formulated as a discrete-time Markov process on the state 
space S = {0, 1, 10, 11, 100, 101, 110, 111, 1000, . . .} where and 1 correspond to unoccupied and occupied sites, 
respectively. In particular, denotes the state that there is no particle on the chain. Note that we do not write 
infinite 0s located left to the leftmost 1. We define |r| by the length of a sequence r. In particular, for each 
element r of S, |t| stands for the length of the system (the position of the leftmost occupied site). Let us write 
the probability of finding a configuration r at time t as P(r;t). The master equation governing the model is 
written as 

P(r; t + 1) = TP(t; i) = £ W(t' -> t)P(t'; t), (6) 

t'gS 

where T is the transition-probability matrix of the process and W(t' — > r) is the transition probability from r' 
to r within one time step. In particular, for r = and r = 1, 

P(0; t + 1) = (1 - a)P(0; t) + (1 - a)0P{l; t), (7) 
P(l; t + 1) = aP(0; t) + (1 - a)(l - /3)P(1; t) + (1 - a)pP(10; t). (8) 

It is difficult to write down the action of T explicitly for the general configuration r because of the parallel 
updat^|. For example, for r = 101101, 

P(101101;i+ 1) 

= (1 - a)(l - p) 2 (l - £)P(101101; t) + (1 - a)(l - p)pP(101110; i) 
+ (1 - - p)(l - /3)P(110101; i) 

+ (1 - a)p(l -p)(l- /3)P(1001101; t) + (1 - a)p 2 P(110110; t) (9) 
+ (1 - a)p 2 P(1001110; t) + (l- a)p 2 {l - p)(l - ^)P(1010101; i) 
+ (1 - a)p 3 P(1010110; t) + ap 2 P(10110; t) 
+ ap(l-p)(l-0)P(lOlOl;t), 

see Fig. El and for r = 100100, 

P(100100; t + 1) 

=(1 - a)(l - p) 2 P(100100; t) + (l- a)(l - p)pP(101000; i) 

+ (1 - a)p(l - p)P(1000100; t) + (1 - a)p 2 P(1001000; t) (10) 
+ (1 - a)(l - p) 2 /3P(100101; t) + (1 - a)(l - p)p i 8P(101001; t) 
+ (1 - a)p(l - p)/3P(1000101; t) + (1 - a)p 2 /3P(1001001; t), 

see Fig. HI 

Note that this model is not equivalent to the parallel-update TASEP with the ordinary open boundary 
condition (open TASEP) [10] . see Fig. [SJ In the open TASEP, particles can enter the system at the fixed 
leftmost site of the finite chain only if it is empty. This means that the open TASEP is a call-loss system. On 
the other hand, in our model, particles can always enter the system. 

4 It is not so difficult in the continuous-time case, see (22) in pQ. 
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Figure 3: List of all configurations which transit to fOf fOf with a non-zero probability within one time step. 



- a n l -p n l ~p 



- a n l ~p n l ~p ^ 



-a Cl l ~P rs p 



-a CI x -p r~* P 



l-a 



n l ~p 



l-a 



l-a f^* 1 



l-a /"X 



rx p r> p 



Figure 4: List of all configurations which transit to f OOf 00 with a non-zero probability within one time step. 



P 



P 



a 





□ 








□ 







Figure 5: TASEP with the ordinary open boundary condition. 



3 Stationary state 

A stationary state (stationary distribution) is a solution to the balance equation 

P(r) = TP(t) (Vt g S) (11) 

with the constraint 

£P(r) = l. (12) 
Let us begin with making an assumption that the stationary-state solution has the following form: 

pm =|, (13) 

1 / a 



where /l is a function on {0, 1} L and satisfies 

/i(lOlr) = /l-i(Ht) +p/ L _ 2 (lr), (15) 

/ i (100T) = (l-p)/ i _i(10T), (16) 

/£(01t)=p/£_i(1t), (17) 

/ l (00t) - 0. (18) 



From P(0) = TP(0), we get 



From P(l) = TP(1), we get 



Thus, 



p(1) 4rt^)|- < 19 > 



pm = 1(^^)1 (20) 



Z \p(l-a)J /T 



/!(!) = |. (21) 
/ 2 (10) = |. (22) 

We reduce the balance equation for the general configuration to that for the open TASEP. One can obtain, for 
all 11t e S with length L, 



/ L (llr) = V -L-W{1W -+ \\T)f L (lW) 
* — ' 1 — a 

T'eAn(r) 

+ V -W(1t' llr)/ L (0lT'), 



(23) 



r'6^lii(r) 



for all 101r G S 1 with length L, 

/i(lOlr) = ]T ^-^(lOlr' 101r)/ L (101r') 

t'£^ioi(t) 

+ V -J— WHO/ -> 101t)/ l (110t') 
1 — a 

t'EAioi(t) 

+ -W(10t' -> 101t)/l(010t / ), 



(24) 
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and for all 100r € S with length L, 

/ L (100r) = ]T 



T'e^lOo(T) 



1 - a 



W(10r' ->■ 100r)/ L (10r') 



E 



i 



(l-a)(l-p) 



W(10t' -> 100t)/ l (00t') 



(25) 



t'G.4i 00 (t) 

see Appendix lAl for detailed calculation. The sets A...(t) are defined as 

„4 n ( T ) :={t'\W(11t' -> llr) > A |r| = |r'|}, (26) 

Awi(t) :={t'\W(W1t' -> 101r) > A |r| = |r'|}, (27) 

^ioo(r) :={t'|W(10t' -> lOOr) > A |r| = |t'| - 1}. (28) 

For example, 

Au(0110) = {0110, 1010,0111, 1011}, (29) 

Aoi(110)={HO,lll}, (30) 

ylioo(llO) = {0110, 1010,0111, 1011}. (31) 

Equations (|2"3")) - (|2~51) give the balance equation for the open TASEP on the L-site chain with the entry probability 
equal to 1 if we regard /l(t) as the probability of finding a configuration r. Thus, if we set to be proportional 
to the stationary-state solution to the open TASEP, the form (fT3|) and (|T4| gives a stationary-state solution to 
our model. Substituting the matrix product solution to the open TASEP which was found in [10] into (|14[) , we 
obtain 



P(r) = 



{W\DX T 



■X T1 \V) t = 1t L - 1 ...t 1 , 



(32) 



where X\ 



v p(l-a) / 

D and X$ = E are matrices, (W\ is a row vector and \V) is a column vector. They should satisfy 
EDEE =(1 - p)EDE + EEE + pEE, 
EDED =EDD + EED + pED, 
DDEE =(1 - p)DDE + (1 - p)DEE + p(l - p)DE, 
DDED =DDD + (1 - p)DED + pDD, 
DDE\V) =(1 - f3)DD\V) + (1 - p)DE\V) + p(l - 0)D\V), 
EDE\V) =(1 - P)ED\V) +EE\V) +pE\V), 

(W\DEE ={l-p) (W\DE, (33) 
(W\DED =(W\DD +p(W\D, 

p(l-/9) 



DD\V) 
ED\V) 



jE\V), 



-D\V), 



(W\EE =0, 
(W\ED =p(W\D, 

which is the algebraic relation found in [10] with a = 1. This algebraic relation is compatible with the assumption 
(p~5|) - (jT5|) . A representation of the relation is given in Appendix [Bl In view of ([2~Tj) , we have to chose the 
normalization as 

P 



(W\D\V) = 



Note that, from the algebraic relation, we can derive (W]D-E|\^) 
normalization constant Z is expressed as 

L 



which is compatible with 



Z = 1 



E 

L>1 



p(l - a) 



(W\D(D + E) L - 1 \V) 



(34) 
The 

(35) 



G 



when the right-hand side converges. This form is similar to the generating function of the normalization constant 
of the open TASEP. The function 



=1 + H L (W\CD((D + E) L ~ l \V) (36) 



L>1 



is useful to calculate the average length of the system and the average number of particles. The case where 
£ = p (-^L a ) an d C = 1 corresponds to the normalization constant: 



z =HiRT^)'V- (37) 

We can easily solve the balance equation of the usual M/M/l queue recursively. (The equation is just a 
three-term recurrence formula.) In our model, however, we have not found a recursive way to solve the balance 
equation. In this sense, our assumption (fT3f and (fT4|) is truly an Ansatz (or working hypothesis). 

In the continuous-time case, a similar form gives its stationary state except that is replaced by a, and an 
overall constant § appears instead that the first D in the matrix product disappears, see (25) in [T]. However, 
the continuous-time limit (first replace a — > a At, f3 — >• j3At and p pAt, and then take the limit At — > 0) of 
the solution (|3"2"|) is exactly the same as the solution to the continuous-time model. 

At this stage, we do not know if the power series for Z (the right-hand side of (j3"5|) ) converges and the solution 
(|3"2"|) is meaningful. Note that the form (j3"2")l without the normalization constant always gives a stationary 
measure of our model, which is a solution to the balance equation (fTTj) without the constraint (fT2"|) . If a process 
on a countable state space is irreducible and recurrent, and has two stationary measures \i and v, then /x is 
proportional to v (Proposition II. 1.3 of [23). Furthermore, if an irreducible process has a stationary state, the 
process is (positive) recurrent and thus its stationary state is unique (Theorem II. 1.1 of [23]). Our process is 
easily seen to be irreducible for generic a, f3 and p. Thus the stationary state (|32p is unique if the power series 
for Z converges and there is no (normalizable) stationary state if it diverges (Corollary II. 1.2 of [23 ). In the 
next section, we obtain the condition on the parameters such that the power series for Z converges and a more 
explicit closed form for Z . 

4 Phase diagram and average values 

We can derive the following closed form for Z(£, £), see Appendix IB1 for details: 

Z( «^T^|ko' (38) 

where 

vff n _ m 1 o ~ e(i - o - v we, o - m - o) 2 - ^(Wxj ,„ Q , 

y K>u- 2(1 +kc) ' ( ^ J) 

*(£,0=(i+pO(i+rtO- (4°) 

The power series for Z(£, 1) (the right-hand side of (|36p with £ = 1) converges to — - s)2 4; when 

2 

Replacing ^ by a ) ; we obtain the region where the power series for Z (the right-hand side of ([35])') converges 



as 



a<a c =^ P<l-^T—p, (42) 
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Table 1: Comparison with the continuous-time model. In the second row, the settings of the models are given. 
In the third row, the stationary currents of the open TASEPs in the limit where the length of the chain L — > oo 
are described. In the fourth row, the critical lines of the exclusive queueing processes are described. The right 
column is obtained by taking the continuous-time limit of the left column. 



discrete time (parallel update) 


continuous time 


entry probability a 


entry rate a 


exit probability P 


exit rate P 


hopping probability p 


hopping rate p 


J do 


J DO 



2 

a(p — a) 
p — a 2 

Pip ~ P) 
p- p 2 



a,0>l- 



a < min(l — y/1 — p, j3) 
P < min(l — y/1 — p, a) 




a,P>% 

- u/p) a < rnin(§,/3) 
P{l-P/p) /3<min(|,a) 



P(p~P) 
p-p 2 



P>1-V^~P 
P<l-V^~p 



P/p) 



P>\ 
P<%. 



and its closed form as 

2(1 -a)p 



R-p + 2(l-a)P' 



(43) 



where R = y/p(p — 4a(l — a)). 

As we see Fig. [6J the region (|42|) is embedded in the region a < P where the usual M/M/l queueing process 
converges. The critical line a — a c consists of two parts; a curve and a straight line. Mathematically, the curve 
is given by a solution to the denominator of the form (1431) = 0, and the straight line by a solution to R = 0. 
Physically, the two parts correspond to the stationary current of the open TASEP in the high-density phase 
(P < min(l — yT — p, a)) and the maximal-current phase (a, P > 1 — y/1 — p), respectively, see Table [T] This 
property is due to the form ([33)1 . In fact, the stationary current of the open TASEP on L-site chain with the 
entry probability 1 is given by 

J L - (44) 

CL + l +PCL 

where Cj, is the coefficient of £ L in Z(£, 1). Note that the critical line for the continuous-time version of the 
model can be also written in terms of the stationary current of the continuous-time open TASEP [T] , see Table 
[T] again. 

Let us calculate some average values in the stationary state. We can calculate the average length of the 
system (the position of the leftmost particle) and the average number of particles in the system by differentiating 
lnZ(£, C) with respect to £ and £, respectively, as 



-p(l-a) 



(N)=^-lnZ{—^X 



= ap{R-p+2(l-a)) 

R{R-p + 2{l-a)PY K ' 

a(l - a)(p - Zap + R) , 

c=1 R(R -p + 2(1 - a)P) ■ [Q> 



The number N q of particles in the queue (i.e., N q — #{j > 2\rj = 1} for t^t^-i ■ ■ -n G S and N q = for 
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Figure 6: Phase diagram for the parallel-update exclusive queueing process. The critical line a = a c (solid 
line) divides the parameter space into the convergent phase and the divergent phase. We drew the dashed 
line representing the critical line of the usual M/M/l queueing process as well. Note that there is a region 
a c < (or <)q < j3 where the M/M/l queueing process converges but our new model diverges. 



0) is also one of the typical values in the queueing theory. We can calculate the average of N q as 

(47) 



W =JW E ? WC^KCD + E) L -\D + E)\V) 



ZdC 

s L>2 



ZdC, 



1 + WTpJo) mo — J~ 



(48) 

,C=i 



= {N)-- p . (49) 
Derivation of the second equality is given in Appendix [Cj 

5 Case p = 1 

In this section we treat the model with p = 1, where the bulk hopping rule is nothing but one of the typical 
deterministic cellular automata, i.e., rule 184. This case was analyzed in [3S]. The matrices and vectors have a 
two-dimensional representation 



°-( J o). £ =(o f ).W-('.V^).|V)-(^), (50) 



see [21] . 

In the case where p = 1, each particle necessarily hops if its right neighbor site is empty. Thus, any 
configurations containing at least one sequence 00 vanish in the stationary stateH In fact, we can see 



It a 



P(1t l _i...00...ti) = - (W\---EE---\V} = 0. (51) 

Z \1 — a / 

One can reduce the stationary-state probability to the following simpler expression: for r = tltl-i---t\ 
(tl = 1), 

!0 if t contains a sequence 00 

Ka)' " 1 ;:?"' ■*>-*• <52) 



5 This means that the process is no longer irreducible when p = 1. But the process on the subset of the state space S consisting 
of and all the configurations without the sequence 00 is irreducible. 
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Figure 7: Phase diagram for p = 1. The critical line a = a c (solid line) divides the parameter space into the 
convergent phase and the divergent phase. We drew the dashed line representing the critical line of the usual 
M/M/l queueing process as well. Note that there is a region a c < a < (3 where the M/M/l queueing process 
converges but our new model diverges. 



Thanks to the two dimensional representation, calculating Z(£, £) is an easy exercise: 

1 



1 



/3(1+«C) 



(3 — a — a/3 

The square root vanishes from Z, and thus the critical line loses the straight line part, see Fig. [TJ 

P 

The average length and the average numbers of particles in the system and the queue are simplified as 

„2 



(L) 



a 



(3 — a — af3 ' 



p — a — a[3 



a 



(3{l3-a-aP) 



(53) 
(54) 

(55) 

(56) 



It is also easy to calculate the probability distributions of L and N by expanding Z(£, 1) and Z{ j^— , £) around 
^ = and ( = 0, respectively: 



Prob[the length of the system = L\ 



1 



1 d L Z(£,l) 



(3 — a/3 — a 



Z\\-a) LI d£ L 
Prob[the number of particles in the system = N] 

i i d N z( T ^,o 



/3{l-a) \/3{l-a) 



(57) 



Z N\ dC, N 

(3 — a/3 — a 
13(1 -a) 

j3 — a/3 — a 



(58) 



/3{1- (3 + a/3)(l ~a)\ /3(1 - a) 



a(l - P + a/3) 



N 



N = 0, 



N > 1. 
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6 Conclusion 



In this paper, we introduced an extension of the discrete-time M/M/l queueing process with excluded-volume 
effect (exclusive queueing process). We took the parallel- update scheme. A stationary-state solution was 
obtained in a slightly arranged matrix product form for the parallel-update open TASEP. We found that the 
critical line which separates the parameter space into regions with and without the stationary state is given 
by the stationary current of the open TASEP in the maximal-current and high-density phases. Particularly, 
we should note that the entry rate cannot be greater than the maximal current of the TASEP if the model 
converges. We calculated the average length of the system (the position of the leftmost particle) and the average 
number of particles in the system. These results recover the stationary state, the critical line and the average 
values of the continuous-time version of the exclusive queueing process [I] , in the limit where the time interval 
At — > 0. In this sense, our discrete-time model is a generalization of the continuous-time model. We calculated 
the number of particles in the queue as well. 

When p = 1, i.e. the bulk hopping is deterministic, the matrices and the vectors constructing the stationary- 
state solution have a two-dimensional representation. The probability distributions of the length of the system 
and the average number of particles are both geometric in this case. 

We leave calculations of the density profile and correlation functions in the stationary state as future works. 
In the TASEP with the ordinary open boundary condition, a domain wall theory explains its phase diagram 
successfully |17j . Investigating how our model converges to the stationary state or diverges with a similar 
argument will be also an interesting study. 
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A Reduction of the balance equation 

In this section, we derive Eqs. (|2U)) - (|2"51) . We calculate the action of T on P under the assumption and 
(fT4"|). Let us use a short hand notation a — a , . For a configuration llr with |llr| = L, the action of T is 



calculated as 



ZTP(Ut) 

Z Y W{11t' llr)P(llr') + Z 



Y W{10W ->■ llr)P(lOlr') 




e-4n(r) 



+ z Y w ( lT ' -> Ht)p(It') 



Y W(11t' -> llr)a L / L (llr') + 



Y W(101t' -> 117> l+ 7 l+ i(101t') 




(59) 



+ W( - lT ' ~> lHa L ~7L-i(lr') 



t'S^ii(t) 





where we used 



W(101t' 



llr) =pW(Ut' -> llr) 



p(l - a) 



W(1t' -> llr) 



(60) 



a 
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and 

/ L+1 (101r') = f L (llr') +pf L -i(lr'), ph-i^W) = /i(Olr'). 

From P(llr) = TP(11t), we get ([55]) ■ 

For a configuration lOlr with |101r| = L, the action of T is calculated as 

ZTP(lQlr) 

=2 ^ FF(101r' -> lOlr)P(lOlr') + £ ^ W(110r' -> IOIt)P(IIOt') 

T'G^lioi(r) T'G^lioi(r) 

+ Z W{1Qt' -> IOIt)P(IOt') + Z Y W(1010r' -> IOIt)P(IOIOt') 

r'e^ioi(r) t'G.4ioi(t) 

+ Z ^ W(1001r' -> lOlr)P(lOOlr') 

t-'G^ioi(t) 

= Y W(101t' -> f01r)a L / L (101r') + ^ W(110r' -> f Of T)a L / L (lf Or') 

t'EAioi(t) t'£A 10 i(t) 

+ Y W(1Qt' -»■ IOItJo^-Vl-iCIOt 7 ) + ^ W(1010r' -> f OfT)a L+1 / L+1 (f Of Or') 

T-'e^ioi(r) r'G^lioi(r) 

+ ^ W(1001r' -> f0fT)a L+1 / L+ i(f00fr') 

r'e^li i(T) 

=a L \ V W(101t'^101t)/ l (1Q1t') + V W(110t' 101r)/ i (H0r') 

\ z — ' 1 — a £ — ' 1 — a 

\T'eAioi(r) t'£A 1(> i(t) 

+ V -VFffOr'^ f0fr)/ L (0f0r') I, 
' a / 

r'G^ioi(r) / 



Moi(t) 

where we used 



W(1001t' -> lOlr) =— ^— VFflOfr' -> fOfr), 
1 -p 

W(1010r' -> fOfr) =pW(110r' -> lOlr) = — - "W (10/ -> fOfr), 

a 



and 



/ L+1 (f0f0r') = / L (ffOr') +p/ L _i(10r'), 
p/i-rClOr') = /i(OfOr'), / i+1 (1001r') = (f - p)/ £ (l Of r'). 

From P(101r) = TP(fOfr), we get (J5JJ). 

For a configuration fOOr with |f00r| = L, the action of T is calculated as 

ZTP(fOOr) 

=Z Y W{10t' ->■ lOOr)P(lQr') + Z Y W(WQt' ->■ fOOr)P(fOOr') 

T'eA W a(r) T'eA 100 (T) 

= Y W(10t' -> fOOT)a L / L (fOT') + Y W(W0t' -> f00r)a L+1 / L+ i(f00r') 

t-'G^Iioo(t) r'eAioo(T) 



r ( V — VFffOr' -> 100r)/ i (10r / ) 
\ ' f — a 

\ T'eA W o(r) 

r'e^ 100 (r) V ^ 
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where we used 



W(100t' -> 100r) 



l-p 



W(10r' lOOr) 



and 



/ L+1 (100r') = (1 -p)/ L (10r'), /i(OOr') = 0. 
From P(100r) = TP(IOOt), we get ||25J|. 



(69) 



B Derivation of Q 

According to [TU], D,E, (W\ and \V) can be reduced as 
D x \ ^ ( E x E 2 



D 



E 



, (W| = ( (W 2 | ), |V) 



r» 2 o ; ' v o o 

where matrices Di,D 2 ,Ei and i?2, and vectors (W\\, (W 2 \, |Vi ) and |Va ) satisfy 



|V1 ) 

1^2 > 



Di^i = (l~p)(D 1 +E 1 +p), 
P(l-j8), 



/3 



E 2 D 2 =p(D 1 + E 1 +p), 
E 2 \V 2 )=p\V 1 ), (W 2 \D 2 = <Wi|p. 



(70) 

(71) 

(72) 

(73) 
(74) 



We can easily show that -D,-E, (W| and |V) satisfy the relations ([531 from the relation (j7Tj) - ((T4"]) . In view of 
the normalization (|34[) . we impose (Wi|Vt) = 1. The following gives a representation for these matrices and 
vectors, which will not be used in calculating Z(£, (): 

(Wi\ = a, 0,0, . . . ), (w 2 \ = v^(o, i, -yr^, (-v/r^p) 2 , • ■ • ), 



/ / p(l-p) 



Di = 



l-p V 1 -p 



l-p V 1 -p 

l-p 



D 9 



l-p 



i) 



/ 



Si 



p(i-p) 



l-p 

V 1 -p 1 -p 



V 1 -p 1 -p 



V 



(75) 



l-p 



\Vi) = 

where we arranged the representation (6.6)-(6.11) in [10 . 



( 1 \ 







\v 2 ) 



(3 nr^ 



1-/3V P 



\Vi), 
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To calculate Z(£, £), we imitate the calculation (|55|) for p = 1 as 

o =i + E ^wcd(cd + s^iv) 



L>1 



--l + K(W\D{l-t(CD + E))- 1 \V) 

'fli OW X iXE 2 \ ( \Vi) 



-l+£C((^x| 0y /^ CjD2X l + eCD 2 XE 2 )\\V 2 ) 

=1 + eC ( (Wi|Di + <W 2 |A! )(X|y 1 )+^£ 2 |F 2 }), 
where X -1 = 1 — £C^i — £-^1 — £ 2 CE 2 D 2 . Using the relation (|T3")) . we can eliminate Z?2 and i?2 in X as 

x- 1 = (i- P 2 eO - £C(i + pO£>i - £(i +KC)£f (77) 

Furthermore, using the relation (f74|) . we can eliminate D 2 ,E 2} (W 2 \ and \V 2 ) in Z(£, £) as 

z&O =i + ^COMA +p)x(i + pt)\v 1 > 

=i + (Wi| {eC(i+p£)£>i+p£C(i+K)}*l^i> 

=i + (w x | {-x- 1 + (l +KC) - £(i + p£C)£i} x\V! > ' ' Nf 



=(1+^^X^x1X1^ 



In the last equality, we used the second relation of (|72|) and (Wi|Vi) = 1. Now we borrow an idea from Section 
4.3 of [3J, where the generating function of the TASEP with a single defect particle on a ring was obtained. Set 

D' = aD x +b, E' = cEx + d (79) 

with a = . rr^.b = 1 — \/u(l — p), c = , 1 and d = 1 — \ — and we get 

= D' + E' (80) 
(Wk |£' = <TO, D'|F 1 > = ('^^a + 6)|F 1 )=:/|F 1 ), (81) 

noting the relations (fTTj) and (|72l) . Set u = , and we get 

x- 1 =(i - P 2 eo - + ko(£i + udo 

=g-h{D> + E!), 



(82) 



where 



fl =(l-^ C ) + «l+rfO^, (83) 



Set ui(l — uS) = and we get 



X = i fl- -(£»' + = - (1 - w-B') -1 (1 - uD'Y 1 . (84) 

9 \ 9 J 9 



Finally, we achieve 



Z(L = l - ±J ^{W 1 \ (1 - uE'Y 1 (1 - wC)" 1 |Vi ) = (85) 
3 5(l-wd)(l-w/) 
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which can be simplified as (|38[) and expanded as 



L>N>1 



A (L - k - 1)\(N - l)l(-p) k 



a LN -(l p) L 2JV g [L - N)\{N - k - l)\k\{N - k)\ 



N-l 

X 

fc=0 



E 



P\ k+1 {L-k-2)\ * (k-£+l)-N\ ; 

/3 J (N-k-l)\(L-N-l)lf^ £\(N-l)\ [ P) 



C Derivation of (N q ) 

In this section, we show 

E Z L (W\CD(CD + E) L - 2 (D + E)\V) = (l + C) - 

The left-hand side is calculated as 

E t L (W\(D((D + E) L - 2 (D + E)\V) 

L>2 
L>1 

~ (! - E t L (W\D((D + E) L - 2 E\V) - £(W\D\V) 

L>2 

= Z(Z,0-l _ _ Q(W\D(1 - £(£0 + E))~ 1 E\V) - § 



Recall the formula 



^-^ D + E ^ = {^n 2 x i4 X cd 2 xe 2 

Using this and the relations (|72|) and (|74|) . we obtain 

(WlDil-aCD + E^ElV) 
= (W 1 \{D 1 +p)X(E 1 +p)\V 1 ) 

= (Wl1 CC(TTk) * WTpJo lVl) 

(Wi\ {(1 +p£0 x- 1 } X {(1 +p0 (l - Ell^M) - x- 1 } |V X ) 
~ e 2 C(i+K)(i+KC) 

Finally, noting (Wi|X|Vl ) = j^B, (Wi\Vi) = 1 and the relation (|73J}, we achieve flggj). 
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